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ABSTRACT 

In the framework of brane-world models lowering the Planck scale to the TeV 
range, it has recently been pointed out that small black holes could be formed 
at particle colliders or by neutrinos interactions in the atmosphere. This article 
aims at reviewing other places and epochs where microscopic black holes could be 
formed : the interstellar medium and the early Universe. The related decay chan- 
nels and the propagation of the emitted particles are studied to conclude that, 
in spite of the large creation rate for such black holes, the amount of produced 
particles do not conflict with experimental data. This shows, from the astronom- 
ical viewpoint, that models with large extra dimensions making the gravity scale 
much lower are compatible with observations. 

Subject headings: PACS Numbers - 04.70.Dy (Quantum aspects of black holes), 
96.40 (Cosmic-rays), 11.25.Wx (String and brane phenomenology), 98.80.-k (Cos- 
mology) 



Introduction 

It has recently been shown that black holes could be formed at future colliders if the 
Planck scale is of order a TeV, as is the case in some extra-dimension scenarios (Dimopoulos & 
G. Landsberg 2001; Giddings & Thomas 2002). This idea has driven a considerable amount of 
interest : see e.g. (Kanti & March-Russell 2002; Casadio & Harms 2002; Chamblin & Nayak 
2002; Cheung 2002; Frolov & Stojkovic 2002; Hossenfelder et al. 2002; Cavaglia 2003; Kanti & 
March- Russell 2003). The same phenomenon could also occur (and could be detected) due to 
ultrahigh energy neutrino interactions in the atmosphere (Feng & Shapere 2002; Anchordoqui 
& Goldberg 2002; Ringwald & Tu 2002; Emparan et al. 2002; Anchordoqui & Goldberg 
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2002; Ahn 2003). Such possibilities open very exciting ways to search for new physics, 
e.g. quantum-gravitational effects through a Gauss-Bonnet term in the action (Barrau et 
al. 2004), Higgs production (Landsberg 2002), evaporation into supersymmetric particles 
(Chamblin et al. 2004) or a D-dimensional cosmological constant (Kanti et al. 2005). This 
article aims at clarifying the situation about the possible consequences of an astrophysical 
and cosmological production of microscopic black holes (/iBH) in the Universe. This is 
especially important as it has been argued that a low Planck scale could already be excluded 
by the available data on galactic cosmic-rays. The first section summarizes the general 
framework of TeV Planck scale models and the basics of cross section computations. Section 
2 is devoted to the investigation of black hole production by the interaction of high energy 
cosmic-rays with the inter-stellar medium (ISM). The subsequent decay of those black holes 
and the propagation of the emitted particles are also studied. Section 3 deals with a few 
other exotic possibilities, including black hole relics and thermal creation of black holes in 
the early Universe. Finally, we conclude that scenarios lowering the Planck scale in the TeV 
era are viable from the astrophysical point of view. 



The "large extra dimensions" scenario (Arkani-Hamed et al. 1998; Antoniadis et al. 
1998; Arkani-Hamed et al. 1999) is a very elegant way to address geometrically the hierarchy 
problem (among others), allowing only the gravity to propagate in the bulk. The Gauss 
law relates the Planck scale of the effective 4D low-energy theory Mp\ with the fundamental 
Planck scale M D through the volume of the compactified dimensions, Vd„ 4 , via: 



It is thus possible to set Mr, ~ TeV without being in contradiction with any currently avail- 
able experimental data. This translates into radii values between a fraction of a millimeter 
and a few Fermi for the compactification radius of the extra dimensions (assumed to be of 
same size and flat, i.e. of toroidal shape). Furthermore, such a small value for the Planck 
energy is not artificial and could be somehow expected so as to minimize the difference be- 
tween the electroweak and Planck scales, as motivated by the construction of this approach. 
In such a scenario, at sub-weak energies, the standard model (SM) fields must be localized 
to a 4-dimensional manifold of weak scale "thickness" in the extra dimensions. As shown 
in (Arkani-Hamed et al. 1998; Antoniadis et al. 1998; Arkani-Hamed et al. 1999), as an ex- 
ample based on a dynamical assumption with D=6, it is possible to build such a SM field 
localization. This is however the non-trivial task of those models. 



1. 



TeV Planck scale framework 
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Another important way for realizing TeV scale gravity arises from properties of warped 
extra-dimensional geometries used in Randall-Sundrum scenarios (Randall 1999). If the warp 
factor is small in the vicinity of the standard model brane, particle masses can take TeV 
values, thereby giving rise to a large hierarchy between the TeV and conventional Planck 
scales (Giddings & Thomas 2002; Giddings & Katz 2001). Strong gravitational effects are 
therefore also expected in high energy scattering processes on the brane. 

In those frameworks, black holes could be formed by any interaction above the Planck 
scale. Two elementary particles with a center-of-mass energy yfs moving in opposite direc- 
tions with an impact parameter less than the horizon radius r + should form a black hole 
of mass M pa s/s with a cross section expected to be of order a ~ -nr\. Those values are 
in fact approximations as the black hole mass will be only a fraction of the center-of-mass 
energy whose exact value depends on the dimensionality of the space-time and on the an- 
gular momentum of the produced black hole. A significant part of the incoming energy is 
lost to gravitational radiation. In particular, lower limits on the inelasticity (ratio of the 
mass of the formed black hole to the incoming parton energy) where derived in (Yoshino & 
Nambu 2002) as a function of the impact parameter. Although those bounds where slightly 
increased in (Yoshino & Rychkov 2005), they still can be considered as reasonable estimates 
of the consequences of the inelasticity and where shown in (Anchordoqui et al. 2002) to dras- 
tically reduce the expected number of produced black holes. To remain conservative (and 
except otherwise stated) the computations performed in this article use a unity inelasticity 
to maximise the emitted flux and strengthen the conclusions. 

To compute the real probability to form black holes in a nuclei collision, it is necessary 
to take into account that only a fraction of the total center-of-mass energy is carried out by 
each parton and to convolve the previous estimate with the parton luminosity (Dimopoulos 
& G. Landsberg 2001). The proton-proton cross section can be written as 
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where f a (x a ) are the parton distribution functions (PDFs) and 
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Those PDFs can be inferred from (Martin et al. 1993) or rescaled from (Dimopoulos & 
G. Landsberg 2001) at various center-of-mass energies by taking advantage of the fact that 
they only depend on y = M BH / y/s after an appropriate change of variables. 

Up to now, the different works related with those phenomena have focused on the pro- 
duction of /iBHs at colliders or through neutrino interactions in the earth atmosphere. The 
first approach is extremely promising as it would allow those fascinating objects to be formed 
and to evaporate within the largest and more sophisticated detectors ever built by human 
beings. Many new physics effects could be discovered, related with the dimensionality of 
space-time, with the radii of extra-dimensions, with new particles, with quantum gravity, 
etc. The second approach is based on the higher center-of-mass energy available when ex- 
tremely high energy neutrinos (allowing for a high signal to noise ratio due do the very small 
standard model cross sections) interact with atmospheric nuclei. 

On the other hand, there are other places in the Universe where one should expect yuBHs 
to be formed if the Planck scale lies in the TeV range. Namely, the ISM where high energy 
cosmic-ray interactions occur above this threshold and the very early Universe where the 
temperature could have been arbitrarily high. The specific energy and time scales associated 
with those particular conditions deserve a dedicated study. 

2. Black hole production by cosmic-ray interactions on the ISM 

The number of black holes produced by cosmic-rays on the interstellar medium per unit 
of time, per unit of volume and per unit of mass can be written as 

dN f°° d(7 

2M pl 

where n is the ISM proton density and <J)(Ecr) is the cosmic-ray differential spectrum at 
energy Eqr- Assuming that the ISM is made of 90% of hydrogen and 10% of helium, n can 
be taken at ~ 1.3 cm~ 3 . As we deal with a diffusion process insensitive to local overdensities 
or underdensities, this average approximation holds very well. For Mo ~ 1 TeV, which is 
the order of magnitude of the lowest Planck scale value compatible with experimental data 
(see references in (Kanti 2004)), the threshold energy for cosmic-rays is 5.5 x 10 14 GeV, 
which lies just below the knee. Measurements from the KASKADE experiment (Roth et al. 
2004) where used to evaluate &(E C r), taking into account the different contributions from 
hydrogen, helium, carbon, silicon and iron. The cross sections have been derived, following 
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(Dimopoulos & G. Landsberg 2001), by the method given in section 2. Table 1 gives the 
resulting differential numbers of yuBHs produced and their initial temperatures as a function 
of their mass for two different Planck scales and two extreme numbers of dimensions (D = 5 
is not included as a very low Planck scale with only one extra dimension is already excluded 
on the ground of the dynamics of the solar system). 

It can be seen that the production rate is small, due to the low value of the flux in 
this PeV-EeV range. It is also very sensitive to the /zBH mass because the cross section for 
p + p — > BH is a fast decreasing function of Mbh- Nevertheless, when the volume of our 
Galaxy ~ 10 60 m 3 is taken in account this leads to ~ 10 19 yuBHs created per second within 
the Milky Way for M D = 1 TeV. The D-dependence is very weak and not included in the 
table as uncertainties clearly dominate at this level of accuracy. 
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11 
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3258 
800 




100 M D 


1Q -55 


11 

6 


2443 
371 



Table 1: Production rate and initial temperature of yuBHs for different numbers of dimensions, 
different Planck scales and different masses. 

To investigate the experimental consequences of those possibly formed black holes, their 
evaporation should be taken into account. As already studied for 4-dimensional primordial 
black holes (Maki et al. 1996; Barrau et al. 2002; Barrau & Ponthieu 2004) antiprotons are 
very promising for this purpose as their abundance is both well known and very small when 
compared to matter cosmic-rays (p/p ~ 10~ 5 ) (Donato et al. 2001). The source term for this 
process can be written as: 
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Jm bh =M d Je cr =J±- i dM BH 



dE p dtdV JMbh =m d je cr =^ 

x n^EcRMMB^Ni^-iEp, M, E CR )dE CR dM BH 

arjp 

where f stands for the number of emitted quanta by a black hole of mass M BB and 
temperature T BH '■ <P ~ M BH /(2T BH ) ~ 2nr + M BH /{D — 3). The relative number of quarks 
or gluons of type % is iVj and the fragmentation function for such a parton of energy T BB 
into an antiproton of kinetic energy Ep is given by dgijdEp. This term also depends on 
Ecr as the produced /xBH is Lorentz boosted with respect to the galactic frame. It has 
been determined using the Lund model PYTHIA (Tjostrand 1994) Monte-Carlo generator, 
between the thresholds E min w "f BH T BH (l ~ Pbh) and E max 7bhT bh (~L + Pbh), as shown 
on the example of Fig. 1 where the important role played by the center-of-mass velocity can 
be easily seen. 



3. Decay and related spectra 



Once this source term is established, the produced antiprotons should be allowed to 
propagate within the Galaxy. For this purpose, a two zone diffusion model described in 
(Donato et al. 2001) & (Maurin et al. 2001) have been used. In this approach, the geometry 
of the Milky- Way is a cylindrical box embedded in a diffusion halo whose extension is still 
subject to large uncertainties. The five parameters used are : K , 5 (describing the diffusion 
coefficient K(E) = K f3R 5 ), the halo half height L, the convective velocity V c and the Alfven 
velocity V a . They have been varied within a given range determined by an exhaustive and 
systematic study of cosmic-ray nuclei data (Maurin et al. 2001) and chosen at their mean 
value. The spectrum is affected by energy losses when antiprotons interact with the galactic 
interstellar matter and by energy gains when reacceleration occurs. These energy changes 
are described by an intricate integro-differential equation (Barrau et al. 2002) where a source 
term q\ er (E) was added, leading to the so-called tertiary component which corresponds to 
inelastic but non-annihilating reactions of p on interstellar matter. 



Figure 2 shows the resulting antiproton spectrum around its maximum. The different 
curves correspond to different astrophysical parameters and to the interstellar and "top of 
atmosphere" fluxes. The inelasticity has also been taken into account (Anchordoqui et al. 
2002) as a simple resclaling as those results are extremely weakly dependant upon the de- 
tails of the shape of the mass spectrum. In any case, they remain extremely far below the 
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background which lies at the level of a few 1CT 2 GeV _1 sr _1 s _1 . But it should be pointed 
out that the spectrum is much harder than the background. The reasons for this are quite 
straightforward. First, although quite high for an exotic process, the cross section for yuBHs 
production is small when compared to SM processes : it is of the order of a few hundreds 
picobarns at the threshold whereas the total pp cross section is not far from 100 mb with an 
antiproton multiplicity close to one for y/s ~ 1 TeV. This explains the overall normalization 
around the maximum. Furthermore, the antiproton flux results from the convolution of the 
primary cosmic-ray spectrum with the fragmentation function. For the secondary component 
(antiprotons due du pp interactions), this fragmentation function is not far from a Dirac dis- 
tribution and the spectrum is expected to have roughly the same slope than the primary one. 
For the /iBH component considered in this article, the fragmentation function is substan- 
tially widened by the boost and becomes the hardest function of the convolution equation. 
The resulting spectrum therefore exhibits the same slope and becomes much harder than 
for the secondary component. The confinement effects due to the galactic magnetic field 
induce a softening of the spectra below the knee but the spectral index difference remains 
the same whatever the energy. Due to this important slope difference, the antiproton flux 
due to /iBHs should dominate over the secondary one around 4 x 10 10 GeV. Although this 
could allow for interesting experimental tests of low Planck scale theories in a far future, 
this remains extremely far above the currently highest energies measured (~ 100 GeV for 
antiprotons). 

Another potential concern about the /iBHs production compatibility with data is related 
to the flux of emitted gamma-rays. They result both from the direct emission and from the 
decay of neutral pions. As the number of partonic degrees of freedom dominates over the 
electromagnetic ones in the standard model, this latter component is clearly dominant. It is 
much higher than the antiproton emission (as pions are much more numerous in the fragmen- 
tation process) but so is the astrophysical background. The resulting gamma-ray spectrum 
has been estimated using the PYTHIA generator to compute the pion distribution and there 
decay. It has then been Lorentz transformed and integrated over the matter distribution of 
the Galaxy. It has also been checked, as in (Barrau 2000), that the extragalactic contribu- 
tion is much lower than the galactic one. The gamma-ray spectrum due to /iBHs remains 
lower than the background up to energies of 10 15 GeV (this latter being extrapolated with 
a power law above the currently measured energy range). Such energies, far above the GZK 
cutoff, are meaningless, especially when considering that the primary cosmic-ray component 
necessary to produce the related /iBHs should have much higher energies. 

The conclusion of this study is that whatever the considered emitted cosmic-rays, /d3Hs 
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that would be currently produced in the Universe if the Planck scale is around a TeV are 
fully compatible with astrophysical experimental data in spite of a production rate around 
10 10 GHz per galaxy. 

4. Thermal production of black holes in the early Universe and relics 

It could also be thought that the possible thermal production of /iBHs in the early 
Universe should conflict with observations. If the temperature of the Universe reached values 
above the D-dimensional Planck scale Mp, /iBHs are indeed expected to have formed through 
the scattering of the thermal radiation. The production and subsequent evaporation of such 
small black holes is irrelevant as long as it takes place before the nucleosynthesis : it just 
contributes to the thermal bath. On the other hand, if black holes mostly decay after the 
nucleosynthesis, the entropy released could modify the relative abundances of light elements 
and contradict the observational data (as it is well known for the gravitino problem in 
inflationary cosmology). The mass of a D-dimensional /iBH with lifetime r is given by : 

M BH = M d (M d t)^ . 

Depending on the number of extra dimensions, the energy required to produce black holes 
that will survive during ~ 100 s (i.e. until the nucleosynthesis) is of the order of 10 20 — 
10 23 GeV, which is above the usual quadri-dimensional Planck scale and clearly meaningless. 
It can therefore be safely considered that /iBHs thermally created in the early Universe are 
not expected to induce any substantial changes in the sensitive primordial 6 Li, 3 He, 4 He 
and D abundances. 

A last point to address is the possible formation of stable relics at the end of the 
evaporation process. Based on several different arguments (e.g. Gauss-Bonnet solutions 
in string gravity (Alexeyev & Barrau 2000) or a renormalization group modification of the 
Schwarzschild metric (Bonanno & Reuter 2000)), it has recently been pointed out that, in 
agreement with the "cosmic censorship" hypothesis, black hole relics are expected to survive 
the Hawking evaporation, even in models with extra-dimensions. In the Planck region, the 
semi-classical approach cannot be safely used anymore and quantum gravitational effects 
must be taken into account. Although no reliable theory is yet established, most attempts 
trying to deal with the endpoint of the evaporation process lead to the conclusion that stable 
relics should be formed. Depending on the details of the underlying model, their mass is 
expected to lie around the Planck mass. The possible contribution of those relics to the 
cold dark matter budget has been extensively studied in the framework of 4-dimensional 
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black holes - see, e.g., (Barrau et al. 2003, 2004) - but it should be investigated whether the 
D-dimensional black holes produced by cosmic-ray interactions on the ISM could lead to a 
substantial amount of remnants. The total mass induced by this process in our Galaxy can 
be computed by integrating the formula given at the beginning of section 2 : 

pOO pt pRgal pVgal /7 AT 

KLs = / / 2nRM rel -———dM BH dtdRdy 

Jm d Jt form Jo Jo dM BH dtdV 

where t is the age of the Universe, tf orm is the galaxy formation epoch, R ga i and y g(d are 
the Milky- Way radius and height and M rei ~ M D is the relics expected mass. This leads to 
a tiny value M*°f ics ks 10 12 g which means that relics dark-matter is negligible and therefore 
remains compatible with observations. 



Conclusion and prospects 

This study shows that the concern related with /zBHs created by astrophysical processes 
(either in the Galaxy or in the primordial Universe) which could conflict with experimental 
data (either through the emitted cosmic-rays or through the possible relics) is irrelevant. 
From the point of view of black hole formation, the earth atmosphere and colliders are 
probably the only places where a TeV Planck scale should have observational consequences. 
Theoretical motivations for extra-dimensional models are therefore not currently contra- 
dicted by astroparticle physics phenomenology. On the other hand, some space is still open 
for new phenomena due to fiBB. evaporation in the Galaxy : as their production cross-section 
has no suppression factor and as their coupling with any quantum during the evaporation 
process is purely gravitational, one can expect an abundant production of rare particles. 
Furthermore, because of the large center-of-mass velocity with respect to the galactic frame, 
those particles are expected to lie in the high energy range with a hard spectrum. Although 
this would make a significant difference with the usual background thermal distributions, a 
clear detection seems quite doubtful and would require a specific investigation. 
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Fig. 1. — Fragmentation function into antiprotons for jets of u quarks emitted by a 200 GeV 
black hole produced by a 10 8 GeV cosmic-ray as a function of the antiproton kinetic energy. 
Upper part : in the fiBB. frame. Lower part : in the galactic frame. 
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Fig. 2. — Antiproton spectra due to the evaporation of /iBHs produced by cosmic-ray inter- 
actions on the ISM. The upper set of curves is an upper limit ignoring inelasticity and the 
lower set of curves takes into account the inelasticity effect as estimated in (Anchordoqui et 
al. 2002). In each case, the interstellar and "top of atmosphere" fluxes are plotted. The dot- 
ted lines are the minimum and maximum spectra (when varying astrophysical parameters) 
and the solid line is the mean spectrum. 



